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ABSTRACT 

Context. X1822-371 is the prototypical accretion disc corona X-ray source, a low-mass X-ray binary viewed at very high inclination, 
thereby allowing the disc structure and extended disc coronal regions to be visible. As the brightest (closest) such source, X1822-371 
is ideal for studying the shape and edge structure of an accretion disc, and comparing with detailed models. 

Aims. We study the structure of the accretion disc in XI 822-371 by modelling the phase-resolved spectra both in optical and X-ray 
regime. 

Methods. We analyse high time resolution optical ESO/VLT spectra of XI 822-371 to study the variability in the emission line 
profiles. In addition, we use data from XMM-Newton space observatory to study phase-resolved as well as high resolution X-ray 
spectra. We apply the Doppler tomography technique to reconstruct a map of the optical emission distribution in the system. We fit 
multi-component models to the X-ray spectra. 

Results. We find that our results from both the optical and X-ray analysis can be explained with a model where the accretion disc 
has a thick rim in the region where the accretion stream impacts the disc. The behaviour of the H/3 line complex implies that some of 
the accreting matter creates an outburst around the accretion stream impact location and that the resulting outflow of matter moves 
both away from the accretion disc and towards the centre of the disc. Such behaviour can be explained by an almost isotropic outflow 
of matter from the accretion stream impact region. The optical emission lines of Hen /14686 and 5411 show double peaked profiles, 
typical for an accretion disc at high inclination. However, their velocities are slower than expected for an accretion disc in a system like 
X1822-371. This, combined with the fact that the Hen emission lines do not get eclipsed during the partial eclipse in the continuum, 
suggests that the line emission does not originate in the orbital plane and is more likely to come from above the accretion disc, for 
example the accretion disc wind. 

Key words. Stars: individual: X1822-371 - Accretion, accretion disks - X-rays: binaries 



' 1 . Introduction the X-ray eclipse by the secondary star has been explained with 

the presence of an accretion disc corona (ADC), an extended 
. XI 822-371 (V691 CrA) is a low-mass X-ray binary (LMXB) cloud of hot ionized gas surrounding the NS and an accretion 
^ ; system, consisting of a neutron star (NS) primary that accretes disc which is big enough so that the secondary cannot totally 
, matter from a Roche-lobe filling low-mass secondary star. The block it from our line f sight. The ADC also explains why the 
X} ■ nature of the primary has been confirmed to be a pulsar with ratio of X - r ay to o ptical luminosity in X 1822-371 is atypically 
H ; the discovery of 0.59 s X-rays pulsations (|Jonker & van der Klis| low (L x /L opt ~ 20 (Grif fiths et all 19781) . usually it is at least one 
, 1200JJ). The orbital period of the system is 5.57 hr and it is seen order of magn itude larger for LMXBs). This is because the ADC 
over a wide range of wavelengths: X-rays ( White et al.lll98ll). is thought to be optically thick in X-rays and thus we cannot see 
UV dPuchnarewicz et al. 1995), optical (|Mason et al.||l980|) and the central source and observe only scattered X-rays, 
infrared (IMason & Cordovalll982l) . The basic concept of emis- 
sion from LMXBs is that the X-rays are produced by the accre- The accretion disc in XI 822-371 is considered to be thick 
tion process and the less energetic radiation through the repro- and have non-axisymmetric vertical structure. Variability of 
cessing of the X-rays in the accretion disc and the secondary star, the light curve outside the eclipse is explained to be due 
The orbital inclinat ion of the system has been determined to be to the extended accret ion disc rim dMason & Cordova! 119821: 
i = 82.5° + 1.5° dHeinz & Now akl l2001l) . The partial nature of iHellier & Masoiil 1 1989b blocking our view of the inner parts 

of the accretion disc. From the X-ray pulse arrival times 

* Based on observations obtained at the European Southern Uonker & van derKlisJ fl2001|) were able to determine that the or- 

Observatory, Chile (Program-ID: 077.D-0603(A)) and with XMM- bit of the system is almost circular and that the mass function 

Newton, an ESA science mission with instruments and contributions of the neutron star is (2.03 + 0.03) X 10~ 2 M Q . ICasares et all 

directly funded by ESA Member States and NASA. (120031) detected Bowen fluorescence from the secondary star 
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and constrained the lower limit of the radial velocity semi- 
amplitude of the secondary to be 300 ± 8 km s . This com- 
bined with the radial velocity semiamplitude of the neutron star 
(calculated from the projecte d semi-major axis of the NS by 
iJonker & van der Klisl d2001h ) constrains the lower limits for 
the component mas ses M 2 > 0.36Mfn and M t > L1 4Mq 
dCasares et alj |2003f). However. iMunoz-Darias et al.1 (2005) take 
the analysis further by estimating the correction for non-uniform 
irradiation of the secondary star, which leads to a lower limit of 
Mi > 1.6M©. The orb ital period of X1822-371 has been ob- 
served to be increa sing dHellier et alJll990tlBavless et al.ll201Ct 
Burd eri et alj |2010). which has been explained by the presence 
of a wind or outflow carrying angular momentum away from the 
system. 

We have organized this paper in the following way. Section 2 
describes our optical observations and data reduction. In Section 
3 we present analysis of the optical data. In Section 4 we 
present X-ray observations, data reduction and analysis. Finally 
in Section 5 we discuss our results and present our conclusions. 

2. Observations and data reduction 

We observed XI 822-371 with the Focal Reducer/Low 
Dispersion Spectrog raph (FORS) 2 instrument 
( App enzeller et al.ll 19981) at the Very Large Telescope (VLT) of 
the European Southern Observatory (ESO) on 15 June 2006. 
The detector of FORS2 consists of two 4096x2048 pixel MIT 
CCDs (15/ffli/pixel) an d they were read out in 2x2 binning 
mode dJehin et al.ll2006l) . 

The data were obtained with grism 1400V (4560-5860 A, 
dispersion 20.8 A mirr 1 ) and a 0.7"slit. The slit was oriented so 
that it included a nearby star ~6.7"south of XI 822-371 for com- 
parison in order to monitor slit losses and possible atmospherical 
changes. The exposure time of a single spectrum was 90 sec- 
onds which, adding the overheads, yielded a time resolution of 
about 2 minutes. In total we obtained 224 spectra covering ap- 
proximately 1.6 orbital cycles of XI 822-371. The seeing varied 
between 0.49"and 1.41"during the night, and averaged 0.87". 

The calibration frames were obtained as part of the standard 
FORS Calibration Plan provided by ESO, these included bias 
frames, screen spectroscopic flat field images and an arc lamp 
spectrum all taken in the morning after the observations. The 
data reduction was conducted by using standard IRAlfl/PyRAlS 
packages for debiasing, flat-fielding, extracting the spectra and 
wavelength calibration. 

Arc lamp spectra were extracted with the same apertures as 
used for the spectra of XI 822-371 and the comparison star. The 
wavelength solution was found by fitting 7 lines with a first order 
cubic spline yielding an RMS of 0.02-0.03 A. The dispersion of 
the data is ~ 0.64 A/pixel and the resolution, measured from the 
arc lamp and the night sky [Oi] lines, is ~ 1.4 A. As the ESO 
calibration plan provides only one arc line exposure, we used 
the position of the sky [Oi] emission line at 5577 A to check 
the wavelength calibration zero point. The line was fitted with 
a Gaussian profile and shifts were calculated as the deviation of 
the line centre from 5577.34 A for both the apertures of X1822- 
371 and the comparison star separately. The shifts were found 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 

2 PyRAF is a product of the Space Telescope Science Institute, which 
is operated by AURA for NASA. 
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Fig. 1. Normalized average spectrum of X1822-371. The rising 
feature on the red end of the spectrum is an artefact from reduc- 
tions at the edge of the CCD. 

to be between -0.465 A and 0.345 A for XI 822-371 and be- 
tween -0.484 A and 0.35 1 A for the comparison star. These shifts 
were then applied to the spectra with the IRAF task SPECSHIFT. 
Naturally this procedure does not take into account the possible 
non-linearity of the wavelength scale and its variation over night 
due to instrumental flexure, but these are considered to be very 
small effects. 

Next the spectra were exported to ASCII format and fur- 
ther processing and analysis were done with the MOLLY soft- 
ware package developed by Tom MarshQln MOLLY the spectra 
were rebinned to the heliocentric frame, then normalised by fit- 
ting a fourth order spline to the continuum. The corresponding 
spline was subtracted from each spectrum. Lastly, the resulting 
continuum-subtracted spectra were binned into 80 ph ase bins us- 
ing the quadratic ephemeris by Baptista et al. (2002). The num- 
ber of spectra per phase bin varied from 1 to 4 (7 bins included 
only one spectrum). 

3. Analysis 

3. 1 . Average spectrum and light curve 

The average spectrum of XI 822-371 is shown in Figure [TJ The 
most prominent feature in the spectrum is the He n 4686 A emis- 
sion line that shows a double-peaked profile typical for an ac- 
cretion disc at high inclination. Another, similar but weaker, 
Hen line is seen at 5411 A. A blend of Cm and Nra lines, the 
Bowen blend, is present in the 4640-50 A region. H/J at 4861 
A has a complex profile composed of superimposed absorption 
and emission components. Other lines that can be seen in the 
average spectrum, which are however difficult or impossible to 
distinguish in individual spectra, we identify as follows. The 
feature at the blue edge of our spectrum could be He n line at 
4541 A detected als o bv lCowlev et all (120031) . but not noted by 
ICasares et al.l d2003l) . The bump on the red side of H/J could be a 
blend by N V at 4945 A and the broad e mission line at 5290 A is 
possibly caused by O vi, as suggested bv lCowlev et"ail d2003l) . 

We also made a light curve by integrating the spectra. The 
values in the light curve are just sums of the counts in individual 
spectra, divided by the corresponding values of the comparison 

3 http://deneb.astro.warwick.ac.uk/phsaap/software/ 
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Fig. 3. Phase binned spectra (10 bins) over the orbital phase. The dotted lines mark the wavelengths of the lines (in case of the 
Bowen blend 4640 A.) 




Orbital phase 



Fig. 2. Folded and binned light curve of X 1 822-37 1 in the wave- 
length range of 4560-5860 A integrated from individual spectra. 



star to correct for the slit losses and variable atmospheric condi- 
tions (this includes an assumption that the comparison star is not 
variable and that it is equally well placed on the slit all the time). 
The errors of the light curve were estimated simply using photon 
noise. Then the light curve was folded and binned into 80 phase 
bins. The result is shown in Figure|2] 



The light curve shows similar structure to previous optical 
light curves. The eclipse is broad and V-shaped and the ingress 
is slightly slower than the steeper egress. The variability outside 
the eclipse is also similar to previously published light curves. 
The slow rise after the eclipse is thought to be due to a thicker 
structure in the accretion disc. However, when interpreting de- 
tails of our light curve one should remember how it has been 
made and take into account the possible sources of error. 

A light curve constructed from the wavelength range of 
5052-5268 A where there are no discernible emission or absorp- 
tion features has similar shape as the light curve in Fig. [2] 



3.2. Spectral line variability 

The variability of the emission lines of the Bowen blend and 
He ii 4686 A, He n 541 1 A and H/3 is shown in Fig.[3]where each 
spectral feature is plotted at 10 different orbital phases. As one 
can see, the Bowen blend is quite faint compared to He n 4686, 
however, there are some changes in its structure. At phases 0. 1 
to 0.4 and 0.6 to 0.9 there are clearly two or three peaks (which 
we identify as 4634.12 A and 4640.64 A from Nm and 4647.4 
A from Cm), with three peaks most prominent at phases 0.6 to 
0.8. The Bowen blend emission is t hought to originate fro m the 
heated face of the secondary star (ICornelisse et al.l l2008). The 
peaks are seen when the secondary is visible to us, that is before 
and after orbital phase 0.5. This implies that around orbital phase 
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0.5 itself there is structure in the disc that is partially blocking 
our view of the donor star. 

The Hen lines at 4686 A and 5411 A show similar vari- 
ability, even though it seems that in the 541 1 A line the central 
trough between the peaks is deeper than in the 4686 A line. More 
analysis of the Bowen blend and He n lines is presented below 
when Doppler tomography is discussed. 

The Hy6 variability is more complex. It appears that there 
is deep phase dependent absorption superimposed on a double 
peaked emission The absorption feature seems to be sweep- 
ing/travelling through the emission from red to blue starting 
around orbital phase = 0.1 - 0.2. The absorption is deepest 
at phase 0.6-0.7, but it is not clear whether this is due to an in- 
crease in the absorption or a decrease in the emission component 
that is "filling" in the profile. 

To study whether the emission lines also get eclipsed during 
the photometric eclipse we integrated the line fluxes (total counts 
above continuum, normalised with the continuum level of the 
comparison star spectrum to correct for atmospheric and other 
changes, errors were estimated again using the photon noise) and 
found that the He n 4686 A and 541 1 A lines do not get eclipsed. 

The integrated He n flux shows clear variability with the or- 
bital phase, see Figure [4] It brightens during the eclipse, which 
can be explained, at least partly, with the decreasing level of the 
continuum. 

The Bowen blend flux does not vary significantly outside the 
eclipse (see Fig.|4]i and actually it resembles the light curve. This 
is well in accordance with the Bowen emission originating from 
the heated face of the secondary star. 

We also checked the behaviour of the spectral feature of O vi 
at 5290A to compare it to the result by Hutchings et al. (2005). 
They find that the O vi emission line at 1032 A shows variability 
over the orbital cycle. Because the feature at 5290A is faint, we 
divied the orbital cycle to only 10 phase bins. The integrated line 
flux shows variability with the orbital phase, but is not eclipsed 
and peaks at phase 0.45, where it is about 1.5 times brighter 
than at other phases. The absence of eclipse suggests that the 
emis sion does not arise close to the compact object, contrary to 
what Hutching s et al.l d2005l) suggest (they have only four data 
points and none at phase 0). 

3.3. Radial velocity curve 

In order to investigate the system parameters of XI 822-371, we 
used the Hen 4686 A line to make a radial velocity curve. This 
was done by fitting two Gaussians to the emission line wings, as 
described bv lSchneider & Yound d!980t) . Because the emission 
line wings are thought to originate from the inner accretion disc, 
their motion will reflect the motion of the neutron star. The width 
of the Gaussians was set to 100 km s _1 and the separation of the 
two Gaussians, a, was varied between 400 and 1500 km s _1 first 
with a step size of 100 km s _1 , later with a step size of 25 km s _1 
around the minimum. Next the equation 



0.5 



ntegrated Hell flux 



V = j - K sin(27T((f> - (p Q )) 



(1) 



was fitted to the resulting radial velocity curves, where V is the 
radial velocity, y is the systemic velocity, K is the radial velocity 
semiamplitude and <p is the orbital phase and </>o the orbital phase 
zero point. 

The results of the fits are presented in a diagnostic diagram in 
Fig El The diagram is used to find the largest usefu l separation of 
the Gaussians as presented bv lShafter et al.l d 1986b . As Figure|5] 
shows, the minimum of cr K /K lies close to a — 1000 km s _1 . We 
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Fig. 4. Integrated line flux of He n 4686 A and 541 1 A {top) and 
the Bowen blend (bottom) normalised with the comparison star. 



elected to use the value of a — 1050 km s based on the lowest 
sum of squares of the fit. With that separation of the Gaussians, 
the resulting fit parameters are: systemic velocity y = —58.5+1 .7 
km s _1 , radial velocity semiamplitude K = 82.9 + 2.5 km s , 
and the orbital phase zero point 0o = 0.06 + 0.004. The radial 
velocity curve (data and the fit) is presented in Fig. [6] The results 
vary little with different values of a, e.g. with a = 950 km s _1 : 
y = -56.1 + 1.6 km s" 1 , K = 76.6 + 2.3 km s"\ and <f> = 
0.76 + 0.004, or with a = 1200 km s _1 :y= -69.2 + 3.0 km s _1 , 
K = 86.9 + 4.3 km s" 1 , and <p = 0.05 ± 0.007. 

These values are in line with previous 
ICasares et all d2003l)_(y _ -43 ± 4 km s _1 , K = 



s" 1 ) and lCowlev et al.l (120031) (y = -106 + 2 km s" 
km s _1 ). 



studies by 
63 + 3 km 
K = 11 + 4 



3.4. Doppler maps 

Although our view of XI 822-371 doe s not fulfil all the axioms 
of Doppler tomography (Marsh 2005), for example the visibil- 
ity of all elements of the disc does not remain constant in the 
case of a structured and thick accretion disc rim at high inclina- 
tion, we applied this t echn ique to our data, as pr eviously done by 
lHarlaftis et aD (1 19971) and lCasares etalJ (120031) . We adopted the 
Dopp ler tomograp hy code where the line flux is allowed to vary 
bv lSteegiisI (120031) . For the systemic velocity we used the value 
obtained from the radial velocity curve fitting (-58 km s _1 ). We 
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Fig. 5. Diagnostic diagram for He n 4686 A. 



produced Doppler tomograms from the Bowen blend (4640 A), 
He n 4686 A and 54 1 1 A, shown in Fig . [7] [8] and [9] respectively. 

The Bowen blend average map (corresponding to the "classi- 
cal" Doppler map) shows a strong bright spot, ring-like structure, 
and a larger, more diffuse and fainter emission region. The bright 
spot falls in the location of the secondary star (V x = 0, V y » 300 
kmjT 1 ), also seen in the Bowen blend map by ICasares et al.1 
(2003). At sufficiently high spectral resolution, the Bowen blend 
has been found to be a power ful diagnostic of the motio n of the 
donor star, as first reported bv lSteeghs & Casaresl (120021) for Sco 
X-l. The ring in the map could be from the accretion disc. There 
seems to be no modulation in the Bowen blend flux as there is 
nothing visible in the modulation amplitude map and the cosine 
and sine modulation maps show o nly variation a t negli gible lev- 
els. In the Bowen blend map by ICasares et alJ (l2003h there is 
only a bright spot at the location of the secondary star, but we see 
also emission spread over larger velocities. This could be due to 
the fact that we used the whole Bowen blen d (centr e d at 4 640 A) 
for reconstructing the map whereas ICasares et al.l (120031) made 
their map with only the 4640 A feature of the blend. The trailed 



spectrogram (top left panel in Fig.|7]shows other components of 
the blend. 

The constant-component Hen maps, 4686 A and 5411 A 
have similar structure, but only 4686 A shows modulation. This 
is at the level of 20% and is located in velocity space at the po- 
sition of the expected accretion stream - accretion disc interac- 
tion/impact zone. Curiously, Hen 5411 A does not show this 
feature. 

It is possible to see in Figures|7j[8]and|9]how the"S-wave" of 
the data has opposite phasing in the Bowen and He n maps. This 
is just as we would expect if the Bowen emission originates on 
the heated face of the donor, and the He n emission comes from 
the accretion disc, i.e. the first traces the motion of the secondary 
and the latter that of the primary. 

The He n emission seems to be peaking at around the veloc- 
ity radius of 300 km s _I and the maximum velocity radius of 
significant emission is around 400 - 600 km s _1 . This is a rela- 
tively low velocity for matter in a neutron star LMXB accretion 
disc. 

Thus we calculated what would be the theoretical observed 
radial velocity from the outer edge of a Keplerian accretion disc. 
To estimate the minimum velocity a disc edge can have, we as- 
sume that the emitting material lies in the orbital plane and that 
its distance from the neutron star is the radius of the LI point 
(a disc larger than that cannot exist). Applying the system pa- 
rameters of XI 822-371 we derive an orbital velocity of ~ 500 
- 600 km s (depending on the value of Mi). If the emission 
comes from the inner disc then the velocity should be larger. 
These velocities are clearly larger than our data shows. Thus this 
implies that the emission cannot originate in the orbital plane of 
the system - a fact also supported by the lack of eclipses in He n 
emission line flux. 



4. X-ray data 

We used archival X-ray data of XI 822-371 obtained with the 
EPIC pn instrument onboard XMM-Newton. The observation 
was done on March 6, 2001 jobs-id 01112301 101) and its dura- 
tion was almost 52 ks covering approximately 2.5 orbital cycles 
of XI 822-371. We produced a light curve and phase-resolved 
spectra from t he Eu ropean Photon Imaging Camera (EPIC pn, 
Striid er et alJ (120011) 1 data and high resolution spectra from 
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Fig. 7. Doppler map of the Bowen blend. Top left: observed 
data, top right: predicted data from the model, middle left: non- 
variable component / constant part of the map, middle right: am- 
plitude of modulation, bottom left: cosine, bottom right: sine. 
The units of data and model are counts over continuum level, 
percentage for the modulation amplitude and its sine and cosine 
components. 



both R eflection Grating Spectrometers (RGS, Iden Herder et al.l 
(I2001la 

The data were processed using the XMM-Newton Science 
Analysis System (SAS). To extract the EPIC pn data we used a 
rectangular aperture centred on the source, since in timing mode 
the E PIC pn data has only one spatial coordinate dStriider et alj 
I2001h and the source appears as a vertical stripe in the image 
made from the event list. 

The extracted X-ray light curve, covering the energy range 
of 0.2-10 keV, is shown in Fig.[l0] The background subtraction 
was omitted since the source is relatively bright and our main 
goal was to study the spectra. 

Next the phase-resolved spectra covering the energy range 
0.3-10 keV were reduced with FLAG=0 and PATTERN=0-4. 
The spectra were background-subtracted with background ex- 
tracted from the data using the same aperture and parameters as 
for the source spectra, but the aperture was placed off the source. 
The spectra were grouped with FTOOLS grppha task so that the 
number of counts in each group of channels was greater than 40. 
We obtained altogether 10 spectra over the orbital period, each 
covering 0. 1 in phase, except the eclipse spectrum which covers 




r t 1 1 i 

— 1 000 1000 -1000 1000 

V x (km/s) V x (km/s) 

Fig. 8. Doppler map of He II 4686 A. Top left: observed data, top 
right: predicted data from the model, middle left: non-variable 
component / constant part of the map, middle right: amplitude of 
modulation, bottom left: cosine, bottom right: sine. The units of 
data and model are counts over continuum level, percentage for 
the modulation amplitude and its sine and cosine components. 



only 0.04 of the phase. The orbital phases of the resulting spectra 
are listed in Table Q] 

The RGS spectra, observed in the high event rate mode, were 
also reduced with SAS, using RGSPROC. The spectra in the en- 
ergy range of 0.3-2.8 keV were grouped with FTOOLS grppha 
task. The spectra contain gaps due to poor calibration and dam- 
age in the detectors. 

The spectra were fitted with XSPEC vll.3 dArnaudlll996h . 
Previous studies of X1822-371 with other X-ray observatories 
such as EXOSAT, BeppoSAX, RXTE and ASCA have shown 
that the spectrum is more complex to model than just a sim- 
ple power law and/or blackbody with an emission line profile 
(usually for the i r on em i ssion line at 6.4 keV). The models of 
Hellier & Mason ( 1989): JParmar et alj d2000h : iHeinz & Nowakl 
(l200lh and llaria et al.l (l200lh were tested for the XMM-Newton 
data but they did not fit adequately. 

Our best-fit model of the XMM-Newton data is composed 
of four emission components and two absorption components. 
These are multicolour black body model for an accretion disc 
(DISKBB), emitting hot diffuse gas model (MEKAL), Gaussian 
line profile (GAUSSIAN), a photon power law (PL), a photo- 
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Fig. 9. Doppler map of He II 541 1 A. Top left: observed data, top 
right: predicted data from the model, middle left: non-variable 
component / constant part of the map, middle right: amplitude of 
modulation, bottom left: cosine, bottom right: sine. The units of 
data and model are counts over continuum level, percentage for 
the modulation amplitude and its sine and cosine components. 

electric absorption model (WABS) and a partial covering frac- 
tion absorption (PCFABS), combined as follows: 

WABS • PCFABS • (DISKBB + MEKAL + GAUSS + PL) 

The physical justification for each component is the following. 
The photo-electric absorption model accounts for the interstellar 
absorption in the line of sight to XI 822-371. The partial cover- 
ing fraction mimics the accretion disc rim of XI 822-371 which 
acts as an optically thick wall preventing a fraction of the emis- 
sion from reaching us. The emission components originate from 
the accretion disc or its immediate vicinity. The disc black body 
is a superposition of multicolour black body elements, follow- 
ing the temperature gradient of an accretion disc. The MEKAL 
model of a thermal plasma models the line emission from the ac- 
cretion disc corona and the bremsstrahlung emanating from the 
accreting matter hitting the neutron star surface. MEKAL was 
not able to fit the strong iron Ka emission line at 6.4 keV, thus, 
a separate Gaussian line profile was required. The non-thermal 
emission is represented with a power law model. In the case of 
XI 822-371 the spectrum includes a non-thermal hard tail which 
is due to photons Comptonized in the accretion disc corona. 

First the model was fitted to the strongest spectrum i.e. the 
spectrum covering phases 0.35-0.45. The Gaussian line profile 




Fig. 10. The folded and binned X-ray (0.2-10 keV) light curve 
of XI 822-371. 



energy and its width were fixed to values of 6.4 and 0.1 keV, 
respectively. Subsequently, after a fit to this first spectrum had 
been found, the values of the normalizations of the disc black 
body, MEKAL and power law were frozen and the model was 
fitted to the rest of the spectra, one spectrum at a time. 

The reduced x 1 values for all the fits are reasonable and lie 
between 1.06 and 1.31. The spectra together with the best-fit 
models at two different orbital phases are plotted in Figure QT| 
The parameter values of the best-fit model are presented in Table 
Q] which shows also the estimated errors corresponding to the 
90% confidence range. 

The covering fraction orbital variation of the partial covering 
model is presented in Figure [12] If the parameter mimics the 
vertical structure of the accretion disc it shows that the height 
of the disc increases at phases prior to the eclipse, that is from 
orbital ph ase 0.6 to eclipse . 

Usin dZombeckl (1 19901) and lJimenez-Garate et ail (120021) we 
identified two emission lines in the RGS 1 first order spectrum. 
The presence of the lines was confirmed by fitting a Gaussian 
line profile at the position of the proposed lines, which are He - 
like Ovu at 0.5685 keV and Lya-like Ovin at 0.654 keV. The 
O vn line is actually a complex composed of a stronger inter- 
combination line at 0.5687 keV and a weaker resonance line at 
0.574 keV. The line components were fit with Gaussian profiles 
and no forbidden component at 0.561 keV was found (see Fig. 
EES. 

As presented by Porquet & Dubau (2000) the line ratios of 
He -like ions can be used as plasma diagnostics. The ionization 
process of the plasma can be determined by calculating the ratio 



G = 



f + i 



where /, i and r stand for the relative intensities of forbidden, 
intercombination and resonance lines, respectively. To calculate 
this ratio we used the equivalent widths from the Gaussian fits 
(see Table 0. This way we obtained the value G = 9.2 which 
suggests that the emitting plasma is purely photoionized (G > 4). 
As we detect no forbidden line, its ratio to the intercombina- 
tion line is zero. This f act combined with Figure 8 and 9 in 
iPorquet & Dub au (2000) infer that the electron density of the 
ac cretion disc corona is h igh, of order n e - 10 12 cm 3 . Figure 1 1 
in Porquet & Dubau ( 2000) shows that the plasma is purely pho- 
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Table 1. Best-fit model parameters and 90% confidence errors for X1822-371 XMM spectra as a function of orbital phase. 



Orbital Phase 


0.05- 
0.15 


0.15- 
0.25 


0.25- 
0.35 


0.35- 
0.45 


0.45- 
0.55 


0.55- 
0.65 


0.65- 
0.75 


0.75- 
0.85 


0.85- 
0.95 


0.98- 
0.02 


N H ,„a bs (xi<n 


0.28^ 


29 + ™ 

-0.02 


3O +0 U2 


27 +u - 05 

"•^'-0.04 


25 +au2 

" -0.02 


26 +U U2 

"•^"-0.02 


0-38^ 


39 * 02 

"-> 7 -0.01 


041+0.02 

"•^ 1 -0.02 


"• J1 -0.05 


N H , P cf al , (xW 22 ) 


Q Q l +0.26 
J - yl -0.27 


t Q9+0.23 
J - y -0.24 


4 70+0.21 
^' -0.20 


4 97+0.28 
^•^'-0.27 


4 38 +0 - 27 


4 97+0.28 

^■'-0.29 


o 45+0.24 

-0.28 


3 6 * 26 

J - u -0.33 


4 44+0.29 


3.26*- 


f 


62+ 02 

"• uz --0.02 


61 +0 01 
"• ul -0.02 


59 +0 02 

" ,J7 -0.02 


59 +003 


60 +o m 

"■""-0.02 


061 +o.oi 

"■ ui -0.02 


73+ 001 


75* 01 
"• /J -o.oi 


74+ 01 
"•' mi.oi 


0.84*;°? 


Tm (keV) 


4 on+o.02 

H -^"-0.02 


4 99+0.02 

-0.02 


5 30 +0 01 


5 41+ - 24 

J - H1 -0.20 


5 iq+0.02 


5 i 9+0.02 
J>J -0.03 


4 04* 01 

^•"^-0.01 


T Q7+0.01 

-'•"-o.oi 


3 90 +0 02 

—0.02 


T 0^0.02 

2 - 8/ -o.oi 


norm diskbb (xl(T 3 ) 


4 54+0.54 

-0.56 


4 54+O.54 


4 54+054 


4 54+0.54 

-0.56 


4 54+0.54 

-0.56 


4 54+0.54 


4 54+0.54 


4 54+0.54 


4 54+0.54 

MJ.56 


4 54+0.54 


kT meka i (keV) 


13 +001 


13 +om 
"• -o.oi 


14 +001 


14 +001 
"•^-o.oi 


13 +001 
"• 1J -o.oi 


13 +om 

"• -0.01 


14* 01 
"- 1 mj.oi 


14* 01 
"- 1 M)01 


14* 01 

"- 1 M101 


12* 02 

"• -0.02 


norm lmkal (xl(T 3 ) 


4 49+4.24 


4 49+4.24 


4 49+4.24 
-t-'-i- -2.37 


4 49+4.24 
H - HZ, -2.37 


4 49+4.24 
-2.37 


4 49+4.24 
H - HZ --2.37 


4 49+4.24 
-2.37 


4 49+4.24 
-2.37 


4 49+4.24 
^•^-2.37 


4 49+4.24 
H - HZ, -2.37 


E G m,.i S (keV) 


6.4 


6.4 


6.4 


6.4 


6.4 


6.4 


6.4 


6.4 


6.4 


6.4 


o- gauss (keV) 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


0.1 


norm gauss (xl(T 5 ) 


1 56 +0 - 60 

1 - JU -0.22 


2 55 +0 ' 52 

iJJ -027 


9 cn+0.52 
^■"-0.33 


O O7+0.42 
2 - 8/ -0.43 


•2 i i +0.26 
-'• Ii -0.65 


9 4Q+0.73 
^•^ -0.22 


1 oq+o.20 

w -0.38 


98 +0 22 

"• ?o -0.30 


1.172* 


0* 16 

"•"-0.18 




2 96 +a2 ° 

^• 7U -0.16 


t i i +0.20 
J - 11 -0.16 


3.30!«;? 2 


2.98*$ 


9 7C+0.19 
Z../J_ 017 


2.80*;?° 


3.64** 


3.61* 2 * 


3.68* 23 


4.6i:°s 


norm p i (xl0~ 3 ) 


3.08^ : | 


3.081°;1 


3.08!°;1 


3.08*1 


3-08^52 


3.08!°;| 


3.08*1 


3-08!«;l 


3.08!o°;l 


3.08*;! 


xl 


1.06 


1.06 


1.08 


1.06 


1.07 


1.10 


1.17 


1.31 


1.29 


1.30 



Table 2. The measured equivalent widths of the forbidden (f), 
intercombination (i) and resonance (r) lines of OVII in the first 
order RGS 1 spectrum. 



Line 


Energy 


EW 




(keV) 


(eV) 


f 


0.561 




i 


0.56875 


14.1 


r 


0.574 


1.53 



toionized (the relative strengths of the lines are correct for purely 
ph otoionized and n e — 10 12 cm~ 3 , similar to our line profile). 

iJi et all (1201 ll) used Chandra data of XI 822-371 and their 
analysis of the O vn line ratio yields similar diagnostics as ours. 

5. Discussion and Conclusions 

In this paper we have presented phase-resolved optical and X-ray 
spectra of XI 822-371. 

The optical and X-ray light curves (Fig.l2land[T0b both show 
the secondary star contributing to the extended partial eclipse. 
The eclipse is narrower in X-rays than in optical: in X-rays it 
lasts approximately 0. 1 in orbital phase whilst in optical it lasts 
at least twice that. This can be understood with the difference in 
the emission region size. The out-of-eclipse variability is thought 
to be due to the varying height of the accretion disc rim. In both 
curves there is a plateau-like section immediately following the 
eclipse egress. This can be interpreted as an accretion disc rim 
that has some vertical structure covering the brighter inner disc 
which is opposite to the stream impact region, thereby producing 
maximum brightness around phase 0.5. 

The semiamplitude of our He n 4686 A radial velocity curve 
is K\ = 82.9 + 2.5 km s which traces the motion of the NS. 



If we interpret the bright spot on the Bowen blend map as a sig- 
nature of the companion star its velocity is approximately Ki = 
300 km s _1 . Thus, we can estimate the mass ratio of the compo- 
nents in the system to be q = K1/K2 = 0.28. With a canonical 
NS mass of 1.4 Mq the companion mass would be 0.39 Mpy 
But if the NS mass is as high as 2.32 Mq (Munoz-Dari as et al.l 
120051) . then the companion mass would be 0.64 Mq. 

The Doppler map of He 11 shows low velocities of the emis- 
sion line peaks implying that the emission does not come from 
the orbital plane of the system (assuming a Keplerian disc). In 
addition, the He 11 flux is actually not eclipsed even if the con- 
tinuum light is. This could be explained with a model where the 
emission does not originate from the accretion disc itself, but 
from matter above the accretion disc, for example a wind, that 
still retain s some of the velocit y fie ld of the disc. Actual ly recent 
papers bv lBavless et alj (1201 Ol) and lBurderi et al.l (1201 Ol) suggest 
the presence of a wind or some kind of outflow from the system. 

The Doppler maps of He 11 4686 A and 541 1 A look similar, 
however, there is no variable component in the 5411 A map. 
This is somewhat peculiar, as it would be expected that these 
lines should originate in roughly the same regions, even though 
they arise from different transitions. The He n emission seems 
to be stronger at negative x-velocities, which corresponds to the 
area where the accretion stream from the secondary star hits the 
accretion disc and gets heated. 

An intrig uing model that can explain our data is that for 
CAL 87 by ISchandl et all (1 19971) . CAL 87 is a close binary 
super-soft source system in the Large Magellanic Cloud and it 
is considered to contain an ADC and a similar thick accretion 
disc as in X1822-371. In this model a "spray" of matter is cre- 
ated in the location where the accretion stream hits th e disc. 
Hydrodynamical simulations bv lArmitage & Liviol (Il998l) result 
in a similar model for the accretion stream impact region in the 
case of inefficient cooling. 
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Fig. 12. The partial covering fraction of the best-fit model over 
the orbital period. 



Fig. 11. The 0.3-10 keV spectrum of X1822-371 at orbital 
phases 0.35-0.45 (top) and 0.75-0.85 (bottom). 



The partial covering fraction of the best-fit X-ray model is 
higher prior to the eclipse than at other orbital phases (excluding 
the eclipse). This implies there is an increased level of absorb- 
ing material in our line of sight, which could be the result of 
the"spray". Also the Hf3 feature in the optical spectra exhibits 
the strongest absorption feature at those phases. 

The complex profile of H/3 could be composed of two 
absorption components, one narrow and one wide, on which 
a double-peaked emission component is superimposed. The 
double-peaked component would resemble the Hen profile, 
originating from matter in a circular motion around the NS, in 
the orbital plane or above it. The narrow absorption is caused 
by the "spraying" matter. The absorption is strongest at orbital 
phases 0.6-0.8 when the spray would be in our line of sight to the 
emitting central parts of the accretion disc. If the spray continues 
moving downstream after the accretion stream impact region, it 
could cause a redshifted feature at phases 0.8-1. This absorption 
component is blueshifted (see Fig. [3}, so it means that matter is 
being "sprayed" towards us, i.e. away from the accretion disc. 
However, in addition to the spray away from the disc, a spray or 
stream over the disc is required to explain the redshifted smaller 
absorption feature after orbital phase 0. 1 . 

The spray model is also supported by other studies, e.g. 
|Ji et al.l d201 ll) observe some X-ray line emission to be located 





irgy (keV) 



Fig. 13. The O vn line at 0.5685 keV, composed of the inter- 
combination line at 0.5687 keV and resonance line at 0.574 keV 
shown together with the fit of the Gaussian line profiles. 



at the stream impact region and Harlaft is et al.l (119971) found Fe n 
/16516 in absorption present only at the orbital phase 0.8. 

In connection with the fact that the Hen emission would 
arise from above the orbital plane, we calculated at what height 
should the emission originate in order to be visible throughout 
the whole orbital cycle. With an inclination of 85° the emission 
should come from a minimum height z = 0.2a, where a is the 
distance between the components. If the inclination was 80° the 
minimum height of z = 0.15a is sufficient so that the emission 
region is seen at all phases. 

Another possibility is that the inclination angle is lower than 
has been modelled previously. It would easily explain the ob- 
served lower velocities and the fact that the He n lines do not get 
eclipsed. However, it would be in contradiction to the eclipses 
seen in the light curves. 
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Our multiwavelength study of the accretion disc in XI 822- 
371 strengthens the case for the existence of non-axisymmetric 
vertical disc structure ("flared disc") in interpreting the observed 
properties. This is the case for both the optical and X-ray data 
and these results agree with t he optical and X-ray resu lts for an- 
other ADC source, UW CrB dHakala et al .1120051120091) . The ac- 
tual nature and cause of this vertical structure is not yet properly 
understood. 
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